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=^*Ground /based  digital  ionosonde  observations  of  the  winter  polar 
cap  FAregion  have  been  used  to  demonstrate  that  the  magnetospher 
ically  induced  ionospheric  convection  can_be  measured  for  the 

bottomside  ionosphere-. . A- number  of^2 4-hour^me_asurement made  at 

Thule,  Greenlana^indicate  that  the  drift  direction'iS' "predomin¬ 
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mcters/second.  measurements,  made  at  Goose  Bay,  Labrador 

•fr  show  a  steady  westward  drift  unTIT~Tocal  magnetic  midnight  and  — 

DD  1473  tOITIONOr  I  NOV  II  IS  OBSOLETE  t  .»»  r  A  C  o  TP  rrr, 


SECURITY  C 


IS  RAOC  (When  Dm*  Holered) 


UNCLASSIFIED 


.  '  7(1.  ABSTRACT  (Continued) 

;  then  a  change  to  an  eastward  drift.  The  ionospheric„driftB  ob- 
|  served  at  Goose  Bay-*are  consistent  with  the  expected  sunward 
.  return  flows  of  the  two-^cell  polar  plasma  convection  pattern. 


The  utility  of  data  from  a  network  of  digital  ionosondes  is  en¬ 
hanced  through  automatic  scaling  of  parameters  needed  for  research 
and  radio  wave  propagation  management.  The  values  of  hmF2  deduced’ 
:  by  real-height  analysis  of  automatically  scaled  Digisonde  iono- 
grams  have  been  compared  with  simple  methods  based  on  routinely  1 
I  scaled  ionospheric  characteristics.  Systematic  discrepancies  were  i 
found  between  the  hmF2  values  obtained  from  the  simple  methods  and 
the  real-height  analysis.  Overestimates  of  15-20  km  were  found 
for  the  night  data  from  five  stations  and  low  solar  activity, 
j  Daytime  discrepancies  are  normally  less,  with  80%  showing  agreement 
t  within  £}I0  km.  Analyses,  similar  to  those  presented  here,  are 
i  useful  /for  defining  the  limits  of  the  simple  methods  and  for 
i  gaining  confidence  in  the  automatically  obtained  data  from  the 
j  expanding,  wor ^d-wide  network  of  digital  ionosondes. 
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FOREWORD 


Within  ■•viral  years,  a  world-wide  network  of 
approximately  40  digital  ionosondee  will  be  operating  to 
further  understanding  of  solar-terrestrial  and  radio  physics. 
The  digital  ionosondes  possess  capabilities  to  provide  data 
bases  that  heretofore  were  limited  to  instruments  requiring  a 
significant  capital  investment  (e.g.  incoherent  scatter  radar) 
or  were  labor  intensive  to  obtain  (e.g.  real-height  analyses). 

With  the  maturation  of  the  Digisonde  hardware,  there 
has  been  associated  analytical  work  to  develop  automatic, 
computer-based  algorithms  that  will  provide  reliable  data 
bases  of  the  aforementioned  types. 

The  University  of  Lowell  Center  for  Atmospheric 
Research,  with  cooperation  and  support  from  the  Air  Force 
Geophysics  Laboratory,  recently  contributed  two  journal  papers 
that  document  the  approaches  and  some  results  of  the 
associated  analytical  work.  The  first  paper  summarizes  the 
technique  and  initial  results  of  using  the  Digisonde  to  obtain 
ionospheric  drifts.  Results  are  presented  from  experiments 
conducted  at  Thule,  Greenland  and  at  Goose  Bay,  Labrador. 
These  results  demonstrate  that  a  systematic  study  of 
ionospheric  drifts  is  now  feasible  with  a  network  of  ground- 
based  digital  ionosondes.  The  second  paper  compares  F  layer 
peak  heights  obtained  from  true  height  analyses  of 
automatically  scaled  ionospheric  parameters  to  those  derived 
from  the  routinely  scaled  ionospheric  characteristics. 
Continued  assessment  of  the  computer-based  analytical  methods 
identify  limits  and  provide  confidence  in  the  data  bases  that 
result  from  the  expanding,  world-wide  network  of  digital 
ionosondes. 
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This  scientific  report  includes  the  recent  journal 
articles  mentioned  above i 

Reinisch,  B.  w. ,  J.  Buchau  and  E.  J.  Weber/  "Digital 
Zonosonde  Observations  of  the  Polar  Cap  F  REgion 
Convection/"  Fhysica  Scripta,  Vol.  36/  pp.  372-377, 
1987. 

McNamara,  L.  F. ,  B.  w.  Reinieoh  and  J.  6.  Tang, 
"Values  of  hmF2  Deduced  from  Automatically  Scaled 
Xonograms,"  Adv.  Space  Res.,  Vol.  7,  No.  6,  pp. 
(6) 53- (6) 36,  1987. 
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Digital  lonoaonde  Observations  of  the  Polar  Cap  P  Region 
Convection* 
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Abstract 

Orourul  him'il  drift  tiS»»rv«non*  of  ilif  turner  polio  enp  f'-rtfm  thmv  ihm 
i hr  oiiiimriorf'lirrieHlI)  induced  lomuphrru  ewivieimn  curt  he  met oured  for 
thr  Mtumvide  mmuphere  A  difiml  lonoromir  with  four  tpumd  mi winy 
■mirnnut  npefdicd  m  Thole,  Greenland  ISO'  CGL)  in  ihr  Duppleodnfi 
modi  A  number  of  2S-ht«ur  mimunminh  indicate  ihm  the  drill  dineiion 
(hint!*  Iinearl)  at  a  function  of  lime  m  accordance  ttuh  dir  predicted 
iiniitunttard  contecimn  paiiem  Thr  drift  vrloctiir*  car>  from  .VSliodOOm  » 
Meaturememt  in « subauroral  station  (Cioou  Hr),  Labrador,  nV  COL) 
mih  Ihr  tame  upaced'anirnnat-Doppler-drifi  technique  thott  r  Jteady  win- 
mo rd  drill  until  loeel  magnetic  midnight  und  a  fast  witch-over  ri  that  time 
lo  rn  ri MfMord  drift  w r  eoneludr  dim  ihr  ob«rved  subauroral  drift#  rrr  (hr 
tunttard  rriurn  flottt  of  Ihr  polrr  plasma  convection,  rnd  Ihr  switch-over 
oecurt  *hrn  ihr  nmion  roimrt  from  ihr  dusk  crll  into  ihr  dr ttn  fill, 

1,  Introduction 

Satellite  observation#  have  e#tabli#hcd  ihe  existence  of  a  two 
cell  convection  pattern  in  the  polar  F  region  drift  when  the 
interplanetary  mugnctic  field  (IMF)  ha#  a  touthward  com' 
ponent  [I,  2),  Thi#  convection  pattern  produce#  an  ami- 
sunward  drift  at  the  hfpheit  latitude#  and  funward  drift  at 
lower  latliutde#  in  the  auroral  zones,  The  plasma  flow  pattern 
is  less  predictable  and  more  Irregular  when  the  IMF  has  a 
northward  component  [3,  4],  Wygant  n  al,  (5]  have  shown 
that  the  response  time  of  the  convection  to  changes  in  the 
IMF  is  two  hours  or  more, 

The  University  of  Lowell  und  the  Air  Force  Geophysics 
Laboratory  have  developed  ground-based  observation 
technique#  that  can  monitor  the  plasma  convection  a#  a 
function  of  time.  We  have  begun  to  measure  the  F*region 
drift  at  three  high  latitude  stations:  ( 1 )  Thule,  Greenland.  86* 
CGL.  (2)  Goose  Bay.  Labrador,  65*  CGL.  and  (3)  Argcntiii. 
Newfoundland.  58*  CGL,  In  thi#  paper  we  give  a  brief 
description  of  the  high-frequency  (HF)  radio  technique  used 
to  measure  the  ionospheric  drift  and  present  #omc  of  the 
results  of  the  Thule  and  Goose  Bay  observations, 

2.  Spaced-antenna  doppler  drift  technique 

Spaced-antcnna  HF  observations  of  Ionospheric  drifts  is  an 
effective  method  to  study  the  dynamics  of  the  ionized  atmos¬ 
phere,  Vertically  transmitted  HF  waves  illuminate  a  large 

*  Thi*  paper  **u*  contributed  in  Ihr  ”*CO*Tr,P  loth  Ouadrennial  Inter- 
natm n/  (tolar-Terrestrial  f’hvucs  Symposium",  Toulnute,  IVHb,  and  tul! 
he  included  in  part  II  of  Ihr  Confirmer  proceedings  ffldnors  H  lluliqsisi. 
()  Hm  und  U  vo i,  ’/Mhm 


urea  of  several  hundred  kilometers  diameter  in  the  F  region, 
an  array  of  untennus  receives  the  signals  reflected  from  the 
ionosphere,  For  the  measurements  described  in  this  paper  we 
used  Diglsondc#  [6.  7],  l,c,.  advanced  digital  loriosondcs.  dial 
operated  alternatively  In  the  ionogram  and  the  drift  modes, 
Actually,  the  lonograms  were  spaced  by  five  minutes  and  the 
time  In  between  was  filled  with  a  number  of  1 6 sec  drift 
measurements,  The  individual  antennas  of  the  receiving  array 
arc  multiplexed  at  the  pulse  repetition  rate  (200  Hz),  The  time 
series  received  at  cuch  of  the  antennas  is  Fourier  transformed 
in  real  time  resulting  in  four  complex  spectra  (In  the  case  of 
four  receiving  antennas)  at  the  end  of  each  measurement 
period,  The  spectral  resolution  was  0,1 25  Hz, 
Cross-correlation  [8]  of  the  complex  spectra  from  the 
spaced  antennas  determines  the  angle  of  arrival  for  each 
spectral  component  containing  significant  signal  energy,  As  a 
result  of  this  analysis  one  can  construct  a  sky  map  showing 
the  location  of  each  reflection  point,  or  source  [9],  specified  by 
a  given  Doppler  frequency,  which  defines  the  radial  velocity 
component  of  the  moving  plasma  for  this  source,  The  Doppler 
frequency  it  Is  given  by 

d,  »  l  vh,  S  *  I,  2, ,,,  (source  Index) 

where  v  Is  the  drift  velocity  and  A  the  wave  vector,  An 
example  of  such  a  sky  map  Is  shown  In  Fig,  I,  The  majority 
of  the  sources  are  located  in  the  south-east  with  a  zenith 
angle  of  about  20*,  The  Doppler  frequencies  vary  from 
■f  1,06  Hz  (labeled  I)  to  -0,81  Hz  (labeled  7),  By  assuming 
that  the  observed  Doppler  shifts  are  the  result  of  a  uniform 
bulk  motion  of  the  reflecting  plasma,  one  can  determine  the 
three-dimensional  velocity  vector  v  which,  in  the  least- 
squares-error  sense,  best  represents  the  Doppler  frequencies 
d,  measured  al  the  source  locations clfled  by  A,, 

It  is  important  to  realize  that  the  14  sources  shown  in  the 
sky  mar  of  Fig,  I  existed  simultaneously  within  a  16  s  time 
window,  Thr  existence  of  the  many  sources  was  first  estab¬ 
lished  by  the  spectral  analysis,  and  then  their  location  was 
found  by  cross-correlating  the  antenm.  signals  In  the  spectral 
domuin,  The  dimensions  of  the  receiving  antenna  array  are 
shown  in  Fig,  2  together  with  the  array  pattern  for  a  10  MHz 
signal,  It  is  evident  that  the  angular  resolution  of  the  array 
would  have  been  much  too  low  to  resolve  the  different  source# 
shown  In  Fig,  I,  If  It  had  not  been  for  lli»  preselection  by  the 
Fourier  analysis, 
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.V  Drift  observations  «l  (loose  Hi) 

In  «n  efTori  10  Minis  (hr  diurnal  variulion  ol  I  hi*  I  region  drift 
uc  calculated  the  uvetuge  drift  velocity  for  every  1 5 mm 
interval  The  direction  in/imulhi  end  I  hr  magnitude  of  I  he 
horizontal  drill  component  ore  plollcd  as  function  of  lime  in 
hy*  t  mid  4  Ibr  two  winter  dim  til  Goose  l#«>  In  both  ciivcv. 
the  drift  is  westwurd  before  local  midnight  and  changes  within 
about  one  hour  to  n  predominantly  eastward  drift.  This 
behavior  o  exactly  whnl  it  expected  m  the  pretence  of  the 
two-cell  polar  convection  pattern  Gootc  Hay  ut  fi?*  Cf.ll.  i* 
under  the  sunward  return  flow  (magnetic  wetii  of  the  dutk 
cell  belorc  local  midnight.  At  midnight,  Goose  liny  rotates 
through  t he  llurung  discontinuity  to  the  dawn  cell  where  the 
sounder  now  tense*  the  sunward  flow  (magnetic  cu»t)  of  the 
dawn  cell 

The  typical  velocity  values  are  around  IOOmvs,  but  at 


itrne*  the  velocity  reaches  values  of  JO Urn  s  or  slightly  above 
(not  shown  in  the  figures)  The  vertical  velocity  component* 
are  indicated  in  the  speed  panels  by  +  (upward)  or  -  (down¬ 
ward)  siynt.  in  yencrul  they  remain  below  50  m  s. 

The  exumplc*  shown  here  were  selected  to  demonsirale  the 
midmyht  velocity  reversal  There  are  many  other  days  where 
the  rcvcrt/il  occurs  but  is  not  us  clearly  defined.  Incases  where 
the  velocities  arc  very  small  our  method  produced  unreliable 
result*  What  needs  to  be  done  in  such  situation*  is  to  increase 
the  Doppler  resolution  by  increasing  the  observation  period 
(integration  time)  from  Ihs  to  J2*  or  more,  in  order  to 
separate  the  diflcrcnl  reflection  points, 

At  times  the  drift  direction  chunyc*  unexpectedly,  os  hi  <i‘ 
AST  on  21  Junuury  I9K2  (lriy,  J),  for  uhoul  three  hours  the 
drift  direction  is  westwurd  instead  of  eastward  us  expected  in 
the  morning  hours,  Such  irregular  changes  could  be  (he  result 
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AM  TN  direction  oi/imulhi  sad  ilu  mnpmiudr  (speed l  nr  thr  hormonal  tanemn  nhservaiinn  *i  ihiee  inundmr  frequencies  give*  drill  si  dirterrni 
tmpnntni  o(  ihtdriH  »rt  thmn  k»  Iimihui  of  am  The  verneal  velncine*  height*  t7W»fm  in  ViTUm)  Solid  curve  carmen*  median  mine* 
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F  -  REGION  DRIFT 

DIGI80NDE  OBSERVATIONS  AT  OOOSE  DAY,  LABRADOR 
20/27  JAN  02  10  TO  05  AST 

hi  4  I  -Kcyi'in  |>ri(i  ,ii  liny,  l.atu^dor  !h  27  J u n  V2  IH  lo  (I.*  AST 


dl  .in  expanding  and  xhrinkmg  size  of  the  convection  pattern. 
VA'Iicn  the  pnticm  expand*  the  untinunwurd  (low  of  the  duwn 
cell  may  Ik  over  Goose  Hay  explaining  the  observed  west¬ 
ward  drill 

4,  Drift  observation*  at  Thule 

In  winter  1 983/ 1 984  AFGL's  Airborne  Ionospheric  Observ¬ 
atory  w.is  deployed  at  Thule,  Grccnlnnd  und  ionoynims  und 
drift  obscrvuiions  were  conducted  with  the  Digisonde  on- 
hoard  ihc  aircrafl  und  four  crossed-loop  antennus  on  the 
runway  nearby.  Results  of  a  number  of  24. hour  runs  arc 
shown  in  Figs  4  und  i  Since  Thule  is  close  lo  the  cenle.  of 
the  polar  cap  we  esped  u  continuous  anlliunward  flow  with 
ii-.  direction  rotating  .W"  m  24  h,  In  magnetic  coordinates  the 
roiaiion  is  approximately  linear  in  time  as  Thule  rotates 
around  the  magnetic  pole.  In  the  examples  shown,  the  general 
trend  of  the  drift  direction  is  approximately  untlxunwurd,  On 
o  December  1910  (Fly.  5)  there  arc  some  large  fluctuation*  In 
the  drift  direction  during  ihc  evening  hours  and  they  coincide 
with  the  occurrence  of  »un-ullgriod  arc*  observed  ori  MOO  A 
with  an  all-sky  imaging  phniomclcr  |9]  I’ig  ft  shows  the 


results  of  four  duys  of  observation  at  the  end  of  January  IUS4 
Again,  ihc  drift  direction  closely  follows  the  antisunssard 
direction  which  is  indicated  in  the  figuie,  The  observed  vel¬ 
ocities  vary  from  about  ,100  to  800 m/s  which  is  three  times 
higher  than  the  typical  values  lit  Goose  liny, 


5,  Conclusion 

The  exploraiory  observations  made  at  a  polar  cup  and  an 
auroral  station  confirm  that  the  ground-based  Digisonde 
Doppler  shift  technique  is  capable  of  monitoring  the  high 
latitude  convection  pattern.  The  average  plasma  velocne' 
vary  from  100  to  100 m/s  ut  the  auroral  station  and  from  inn 
to  800  m  s  m  the  polar  cap  Since  the  Danish  Mcicoiologu.il 
Institute  is  now  operating  a  Digisonde  at  Qunaq  (H?*1  CGI.) 
and  the  Air  Force  Geophysics  Laboratory  operates  Digisonde* 
at  Goose  Hay  (6V  CGL)  und  Argcniiu  (.IN'1  CO  I. )  it  will  Ik- 
possible  to  more  closely  monitor  the  polar  convection  and  to 
correlate  the  observations  with  other  geophysical  phenomena 
J 10.  II),  especially  the  geomagnetic  activity,  and  with  the 
variations  In  the  IMF 

r  M>,  S  ufi,  x 

I  > 


Digital  limoxtindf  Observations  of  the  I'nlar  Cup  I  Region  Convection  .377 
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Fif  6  F-Region  Drifi  at  Thule.  30  Jan  to  4  Feb  84.  Direction  of  the  drift 
for  four  winter  days  is  shown  with  reference  to  the  antisunward  direction. 
The  velocity  panel  shows  the  mean  velocity  and  the  variations  during  the 
four  days. 


Acknowledgement 

The  authors  warn  to  thank  Mr.  Claude  G.  Dozois  of  ULCAR  for  the 
processing  of  the  drift  data 

References 

1  CaufTman.  D.  P.  and  Gurnett,  D  A.  Space  Sci.  Res  1.3,  369  (1973) 

2  Rein.  P  H  .  Spiro.  R  W  and  Hill.  T.  W  J.  Geophys  Res.  86.  7639 
(1981). 

3  Rc/hcnos.  Fl  V„  Planet.  Space  Phys  29.  687  (1981) 

4  Burke.  W,  1 .  Kelley.  M  C..  Sagalvn.  R  C.  Smiddcv.  M.  and  Lai. 
S  T..  Geophys.  Res’  Lett.  6.  21  (1979). 


5.  Wygant.  J.  R..  Torbert.  R.  B.  and  Mozer.  F.  S..  J  Geophys  Res  85. 
5727  (1983). 

6  Bibl.  K.  and  Rcimsch.  B  W.  Radio  Sci  13.  519  (1978) 

7.  Rcimsch.  B.  W„  Radio  Sci.  21.  3.  331  (1986). 

8.  Bibl.  K..  Pfisler.  W..  Rcimsch.  B  W.  and  Sales.  G.  S  .  COSPAR. 
Space  Research  XV.  pp  405-411  (1975) 

9.  Dozois.  C.  G  .  AFGL-TR-83-0202.  ULRF-424  CAR.  NTIS  -'ADA 
140509  (1983). 

10  Weber.  E.  J..  Buchau.  J..  Moore.  J.  G..  Sharber.  J  R..  Livingston. 
R.  C..  Winningham.  J.  D..  and  Rcimsch.  B  W,,  J.  Geophys  Res.  89. 
1683  (19841. 

11  Buchau.  J..  Rcinisch  B.  W..  Weber.  E  J  and  Moore.  J  G.  Radio 
Sci  18.995  (1983). 


8 


f‘h\  Wi  <1  Si  tif'1%1  .if' 


•Wi  r  /io  N  o!  **,  \'i»  ft,  pp  (MM  I *#S 7 

1‘mmu^I  in  Him.iih  All  m*ht\  icmmuM 


<»:?>  !  I  IVHl  Vo  00  ♦  ..NO 
C'opvritrht  (£/  COM'AK 


'  "V 


Permission  for  further  roDmrt^l^  L°  roprortuCtt  and  sol!  !!»!*  report. 
thfl  copyright  owner.  °  °n  by  oth*r‘  mu>t  otoftined  from 


VALUES  OF  hmF2  DEDUCED  FROM 
AUTOMATICALLY  SCALED  IONOGRAMS 


Leo  F.  McNamara/  Bodo  \\’.  Reinisch  and  Jane  S.  Tang 


Lhtiversitv  of  Lowell,  Center  for  Atmospheric  Research,  -150  Aiken  Street,  Lowell, 
MA  01854.  U.S.A. 

Wow  with  Andrew  Antennas,  innovation  House,  Technology  Park  5005, 
Australia 


ABSTRACT 

The  values  of  hmF2  deduced  by  real-height  analysis  of  automatically  scaled  vertical  inci¬ 
dence  Dlgisonde  ionograms  have  been  compared  with  those  obtained  using  simple  methods  based 
on  routinely  scaled  ionospheric  characteristics.  Comparisons  have  been  made  for  ionograms 
recorded  at  five  stations  during  low  solar  activity.  The  real-height  analyses  used 
Titheridge's  program  POLAN,  while  the  simple  methods  used  were  those  of  Dudeney  and  Bilitza. 
Systematic  discrepancies  beeween  the  simple-model  and  POLAN  values  of  hmF2  were  found,  with 
15-2C  km  overestimates  of  hmF2  by  the  simple  methods  during  the  night.  The  majority  of 
discrepancies  lie  between  ±20  km.  During  the  day,  802  of  discrepancies  lie  between  ±10  km. 
Significantly  larger  discrepancies  can  usually  be  attributed  to  special  ionospheric  condi¬ 
tions  or  ionograms,  failure  of  the  automatic  scaling  technique,  or  failure  of  one  of  the 
basic  assumptions  of  the  simple  models. 

SIMPLE  MODELS  OF  THE  IONOSPHERE 


Simple  models  of  the  ionosphere  have  been  used  for  many  years  in  applications  in  which  the 
errors  in  the  resulting  calculations  are  accepted  because  of  the  enormous  decrease  in  com¬ 
putation  times  which  such  models  permit.  One  of  the  most  important  practical  applications 
i6  in  the  calculation  of  sky-wave  field  strengths  and  transmission  losses  at  HF.  See,  for 
example,  the  Supplement  to  Report  252  / 1 / .  For  these  calculations,  the  propagation  modes 
involving  reflection  from  the  E  and  F  layers  are  determined  using  an  ionospheric  model  with 
parameters  which  depend  on  the  routinely  scaled  values  of  foE,  foF2,  M(3000)F2  and  h'F,F2. 


The  model  currently  recommended  by  CCIR  consists  of: 

a  parabolic  E  layer  below  its  height  of  maximum  electron  density,  hmE,  with 
semi-thickness  ymE.  hmE  and  ymE  are  set  to  110  km  and  20  km,  respectively. 

a  parabolic  F2  layer  with  height  of  maximum  density  hmF2,  and  semi-thickness 
yraF2. 

a  linear  increase  of  electron  density  between  hmE  and  the  point  on  the  para¬ 
bolic  F2  layer  where  the  plasma  frequency  is  1.7  foE. 

The  model  value  of  hmF2  is  given  by  the  empirical  equations  /1 ,  2/: 

hmF2  -  1490/[M(3000)F2  +  dM]  -  176 

•  !fb  dM  -  0. 18/ (X  -  1.4)  +  0.096  (R12  -  25) /I 50 

and  X  “  foF2/foE,  or  1.7, 


whichever  is  the  larger.  R12  is  the  12-month  smoothed  sunspot  number. 

The  term  dM  is  an  empirical  correction  term  which  takes  into  account  the  effects  of  under¬ 
lying  ionization  not  allowed  for  in  the  original  Shimazaki  /3/  forraulatior.  Model  estimates 
of  hmF2  are  usually  correct  to  within  20  to  30  km  !\ /.  Dudeney  /4/  has  described  an  im¬ 
proved  model  which  uses  a  cosine  F2  layer  shape  and  a  more  realistic  E-F  transition  region, 
and  has  deduced  that  this  model  yields  more  accurate  values  of  hmF2  than  the  Bradley-Dudoney 
model.  Dudeney  gives  a  revised  formula  for  hmF2: 

hmF2  -  1490  F/(M  +  dM)  -  176 
where  dM  -  0.253/(X  -  1.215)  -  0.012 


(h|?3 
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and  F  -  M  [(0.0196  M2  +  0/(1 . 2967  M2  -  l)]0,5. 

The  value  of  X  is  contralned  to  exceed  1.215. 

In  general,  we  have  concentrated  on  the  Dudeney  /4/  model,  but  we  have  also  studied  the 
relative  validity  of  the  Bllitza  et  al.  /5/  and  Bradley  and  Dudeney  111  models.  We  have 
analyzed  about  1000  Diglsonde  256  ionograms  from  five  North  American  stations:  Lowell, 
Massachusetts;  Argentla,  Newfoundland;  Goose  Bay,  Labrador;  Richfield,  Utah;  and  Erie, 
Colorado.  Tho  ionograms  were  chosen  on  the  basis  of  their  ready  availability,  and  do  not 
necessarily  represent  a  statistically  perfect  sample. 

The  required  values  of  foE,  foF2  and  M(3000)F2  were  provided  by  the  autoscaling  software 
ARTIST  which  is  part  of  the  Diglsonde  / 6,  7/,  The  virtual  height  traces  scaled  from  the 
ionograms  using  ARTIST  were  passed  to  the  true-height  analysis  program  P0LAN  /8/,  to  pro¬ 
vide  reliable  estimates  of  hmF2. 

DISCREPANCIES  BETWEEN  POLAN  AND  SIMPLE  MODEL  VALUES  OF  hmF2 

We  have  analyzed  the  differences  between  the  POLAN  and  Dudeney  values  of  hmF2  according  to 
the  value  of  X  (•  foF2/foE)  and  to  the  value  of  the  correction  term  dM.  The  differences 
have  been  distributed  among  bins  5  ten  wide,  with  extreme  negative/positive  values  being 
placed  in  the  <-30/>30  km  bin.  Negative  discrepancies  indicate  that  the  POLAN  value  was 
less  than  che  Dudeney  value. 

Tables  1  and  2  show  the  distributions  of  Che  discrepancies  for  dM  values  of  0.0  to  0.1,  and 
0.1)  to  0.5,  respectively.  Broadly  speaking,  the  low  values  of  dM  (Table  1)  correspond  to 
nighttime  ionograms,  and  the  high  values  to  daytime  ionograms.  It  can  be  seen  that  the 
Dudeney  result  is  systematically  higher  than  the  POLAN  result  at  night,  by  about  15  t  15 
km.  During  the  day,  there  is  no  significant  systematic  discrepancy  between  the  two  sets  of 
results,  and  81%  of  the  discrepancies  lie  within  110  km.  [Note  that  the  extreme  values  are 
ignored  when  calculating  the  percentages.]  There  were  69  Argentia  and  Erie  ionograms  for 
which  dM  exceeded  0.5.  When  the  discrepancies  greater  than  30  km  are  ignored,  the  average 
discrepancy  for  the  remaining  48  cases  was  15  -  20  km,  the  Dudeney  values  being  the  larger. 
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Figure  1  shows  the  distributions  of  the  discrepancies  for  four  ranges  of  X.  The  discrepan¬ 
cies  change  from  being  systematically  positive  for  low  X  values,  to  being  systematically 
negative  for  large  X  values. 


o  -  DUDENEY 


X>  6 


4<X*£6 


2<X^4 


0<X^2 


Figure  1.  Distribution  of  the  discrepancies  in  kn  between  the 
POLA.1;  and  model  values  of  the  height  of  the  Pl-laver 
peak,  hjnFi ,  The  ordir.ites  are  percentages,  and 
values  lying  outside  t.-.e  30  ion  limits  have  been 
ignored.  The  parameter  X  is  defined  bv  X  =  foT’/foD. 
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The  systematically  high  values  of  lmF2  given  by  the  simple  formulas  at  night  cun  be  trucud 
to  the  models'  neglect  of  underlying  ionization.  Thr  effect  of  tills  underlying  Ionisation 
was  originally  considered  to  be  Insignificant  a  decade  or  so  ago,  which  It  Indeed  was  In 
comparison  with  the  effect  of  the  E  and  FI  layers  during  the  day.  However,  the  daytime 
models  have  been  Improved  to  the  point  where  the  nighttime  errors  are  now  the  largest 
obtained. 
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